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Abstract. Choose n random, independent points in R d according to a fixed distribution. 
The convex hull of these points is a random polytope. In some cases, central limit theorems 
have been proven for the components of /-vectors of random polytopes constructed in this 
and similar ways. In this paper, we provide numerical evidence that the components of 
the /-vectors of random polytopes generated according to five different distributions are 
approximately jointly Gaussian for large n. 

1. Introduction 

Phenomena in high dimensions have been the subject of many studies. In particular, one 
may be interested in studying higher dimensional geometric objects called polytopes. A 
polytope P G R d is a bounded subset of R d formed by taking the intersection of finitely many 
halfspaces. The /-vector of a polytope P is a vector / = (/ , f%..., fd-i) such that /; is the 
number of z-dimensional faces of P. In particular, /o is the number of vertices of P and fd-i 
is the number of (d — l)-dimensional faces, known as facets. 

There are some known results on the distribution of the /-vector of types of random 
polytopes. Previous work has found that individual components of the /-vectors are ap- 
proximately Gaussian when the polytopes are generated according to various distributions. 
In particular, Barany and Vu [3] proved that if a polytope is taken as the convex hull of 
independent and identically distributed (i.i.d.) points according to the standard Gaussian 
distribution, then the volume and individual components of the /-vector of that polytope 
satisfy a central limit theorem; that is, converge to a Gaussian distribution as the number 
of i.i.d. random points tends to infinity. Barany and Reitzner [2] proved central limit theo- 
rems for the volume and components of the /-vector of a so-called Poisson random polytope; 
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that is, the convex hull of the intersection of an arbitrary volume 1 convex body K with a 
Poisson process x- Following earlier work of Reitzner [I] , Vu [5] proved that the volume and 
components of the /-vector of a random polytope drawn from points inside a smooth convex 
set also satisfy central limit theorems. 

In this paper, we consider the joint distributions of the /-vectors of random polytopes. 
The method used in this paper for constructing a random polytope is taking the convex hull 
of a set of i.i.d. random points generated from a fixed distribution. Our conjecture is the 
following: 

1.1. Conjecture. Let N be a positive integer. Let Xn = {Xi, X n } be a collection of 
random points Xi G R d , which are independent and identically distributed according to a 
fixed distribution. Letting Pjy £ ^ d be the convex hull o/Xjv, under mild assumptions on the 
distribution of the Xi, the joint distribution of the f -vector satisifes a central limit theorem. 

1.2. Results. In the five underlying distributions of the random points considered, numer- 
ical evidence supports the conjecture. The metric used to compare the closeness of the 
distribution of the /-vectors with the Gaussian distribution is based on the Kolmogorov dis- 
tance, Dk- For now one may think of as a tool that measures the distance between two 
distributions and that < < 1. If is a small value then the two distributions are con- 
sidered to be close to each other. The values of Dx presented in the table are for comparison 
to the appropriate Gaussian distribution; details are discussed in section 3. Here, d is the di- 
mension, n is the number of random points, and iV is the sample size or number of /-vectors. 



TABLE 1. Uniform Distribution in Cube 



d n N 



D K 



5 64000 25000 

6 64000 25000 

7 8000 25000 

8 1000 3125 



0.008995 
0.007590 
0.007133 
0.02613 



2 



Table 2. Uniform Distribution in £ 1 Ball 



d 


n 


N 


D K 


5 


64000 


25000 


0.009971 


6 


64000 


25000 


0.01075 


7 


4000 


25000 


0.01113 


8 


1000 


2000 


0.03003 



Table 3. Uniform Distribution in £ 2 Ball 



d 


n 


N 


D K 


5 


64000 


25000 


0.008495 


6 


64000 


25000 


0.007873 


7 


2000 


25000 


0.01182 


8 


1000 


940 


0.04573 



Table 4. Standard Gaussian Distribution 



d 


n 


N 


D K 


5 


64000 


25000 


0.01197 


6 


64000 


25000 


0.01125 


7 


64000 


25000 


0.009252 


8 


32000 


1000 


0.01530 



Table 5. Uniform Distribution in Hemisphere 



d 


n 


N 


D K 


5 


64000 


25000 


0.006166 


6 


4000 


25000 


0.008827 


7 


500 


12500 


0.02883 



Since Dk is a metric that takes on values in the interval [0,1], a value with order of magnitude 
10~ 3 as the Kolmogorov distance should be considered a rather small value. 

Figure 1 contains histograms of individual components of 64,000 /-vectors computed from 
the convex hull of 25,000 i.i.d. random points uniformly distributed in the 5-dimensional 
cube. After standardizing the /-vectors to have sample mean zero and identity sample 
covariance, the resulting data lie in R 2 . As further illustration, the 2-dimensional data of the 
standardized /-vector is plotted in R 3 . 
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Figure 1. Histograms of components of /-vectors in R 5 . Data came from 
25,000 /-vectors computed from random polytopes generated by 64,000 i.i.d. 
uniformly distributed points in the 5-dimensional cube. 




Figure 2. Standardized /-vectors in R 5 . Data came from 25,000 /-vectors 
computed from random polytopes generated by 64,000 i.i.d. uniformly dis- 
tributed points in the 5-dimensional cube. 



2. Background 

It was mentioned in the introduction that one way to construct a random polytope is 
to take the convex hull of random points. One way to generate random points is to let 

4 



{Xi}™ =1 be a collection of points generated independently according to a fixed distribution. 
Five distributions used to generate random points were considered and the random polytope 
studied was the convex hull of the {Xj}" =1 . The underlying distributions considered are the 
following: 



(1) Uniform in the Cube. Define B^ := {x G R d : \ \x\ |oo < 1} where | |x||oo := max \xi\. A 

l<i<d 

random point X G B^ is constructed by picking each component of X independently 
according to the uniform distribution in the interval [0,1]. The resulting random point 
X is uniformly distributed in B^. 

(2) Uniform in the Euclidean Ball. Define B d in the £2 norm as B d := {x G R d : ||x|| 2 < 

( d \ 3 

1} where ||x||2 := Z \ x i\ 2 ■ To generate a random point X G R d such that 
\i=i J 

X is uniformly distributed in B d , first generate a point Z = (Z 1 , Z d+2 ) that is 
distributed according to the standard d + 2 dimensional Gaussian distribution. This 
is done by picking each component of Z independently according to the standard 

(Z Z) f d+2 t 19 V 

univariate normal distribution. Let X = \\z\\ 2 ' wnere ll^lh = I Z) J • Then 
it is classical that X is uniformly distributed in B d . See [1] for a proof. 

(3) Uniform in the t\ ball. Define Bf in the l\ norm as Bf := {x G R d : < 1} where 
H^lli : = J2 Let Z = (Zi, Zd+i) where the components of Z are generated 

i=l 

independently according to the distribution with the density function f(x) = |e _ ' a '. 
Let X = ; where ||^||i = Z it is proven in [1J that X is uniformly 

distributed in Bf. 

(4) The Standard Normal distribution. Let X G R d be a random vector and generate 
each components of X independently according to the standard normal distribution. 
Then X G R d is distributed standard normal random vector. 

(5) Uniform in the Hemisphere. Similar to (2), but take the absolute value of the first 
component of the random point X in (2). Note that unlike the underlying convex 
bodies in (l)-(3), this body is neither smooth nor a polytope. 
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3. Computation and Analysis of /-vectors 

The software MATLAB and QHULL were used to carry out the following simulations and 
computations. One realization of the /-vector is computed by first fixing the dimension d, 
generating n i.i.d. points {X l }™ =1 from one of the fixed probability distributions described 
above, taking the convex hull of these points to construct a random polytope P, and com- 
puting the /-vector of the random polytope. Assuming the random polytope is simplicial, 
that is, each of its facets has exactly d vertices, the /-vector of the polytope can be directly 
computed from the facets. Let ei, ...e m C R d be the m facets of P, where e it j is the jth vertex 
of the zth facet, for 1 < % < m, 1 < j < d. By the definition of /, fa-i = m. Since /^_ 2 is 
the number of (d — 2)-dimensional faces of P, fa~2 is found by counting all intersections of 
size (d — 1) between pairs of facets and ej, for i ^ j. Since each (d — 2)-dimensional face 
is the intersection of exactly two facets, this accounts for all (d — 2)-dimensional faces. For 
example, 3-dimensional faces (xi,x 2 ,x 3 ) and (xi, x 2 , X4), where each Xi is a vertex for that 
face, have intersection (xi,x 2 ) which is of size 2 so it must be an edge, but (xi,x 2 ,x 3 ) and 
(xi, £4, £5) have intersection (xi) which is of size 1 so it is not an edge. Let g±, ...,gk G R"' -1 
be these distinct intersections of size (d — 2); clearly /^_ 2 = k. By keeping track of all the 
intersections of size (d— 1), we obtain a list of all (d — 2)-dimensional faces of P. We can then 
repeat this process with those faces. Continuing inductively in this manner, the remaining 
components of / can be determined. 

Note that the convex hull of points in general can be non-simplicial (a facet may have 
more than d vertices), and so the above algorithm for constructing the /-vector from the 
facets is not valid for all polytopes. However, the probability of a non-simplicial convex hull 
arising from the distributions of points considered in this paper is 0. 

A collection of iV realizations of /-vectors was obtained, and the data were then stan- 
dardized and compared to the standard Gaussian distribution. The following demonstrates 
the computations of standardizing the /-vectors. Let f k be the k th /-vector; then ff is 
the i th component of the k th /-vector. Define the sample mean / of the /-vectors as / = 
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(/o) hi •••■> fd-i) where fj = ^J2 fj- Given the sample data {f k }^ =1 , the sample covariance 



fc=i 



N 



matrix S = [sij]fj =1 is a d-by-d matrix with entries given by = -^-j- (ft — fi)(fj ~ fj)- 

k=l 

To standardize to identity sample covariance, first diagonalize the sample covariance S = 
UDU T , where D is a diagonal matrix of eigenvalues and U is an orthogonal matrix where 
the columns of U are the eigenvectors. Moreover, the eigenvalues of D can be assumed to 
be in decreasing order, Ai > A 2 > .--Xd- Note that the components of the /-vectors have 
some dependence on each other. For polygons, the components of the /-vector obey /o = /i- 
For convex polyhedra in R 3 , the components of the /-vector satisfy /o — f\ + fi — 2, known 
as Euler's relation. In higher dimensions, linear dependence between the components of the 
/-vector causes the covariance matrix of the /-vector to be singular. These type of linear de- 
pendencies are known as Dehn-Sommerville equations for simplicial polytopes [6]. A singular 
covariance matrix leads to zero eigenvalues in the diagonalization of the covariance matrix. 
In practice, the distinction between non-zero eigenvalues and approximately zero eigenvalues 
of the sample covariance matrix was easy to identify. In such case, let the eigenvalues of 

the sample covariance matrix that are close to zero be eliminated to ensure that is 

Ai ••• 



well-defined. So define a new matrix D* 







A r , 



where D* is a p-bj-p matrix with 



the eigenvalues of matrix D that are close to zero eliminated. Also define a new matrix 



U* 



Ul 



where the eigenvector U; L corresponds to the \ eigenvalue. Then let 



/* = (D*)-^(U*) T f k , where /* is the k th re-centered /-vector. Then {/ fc }f =1 C R p is now a 
standardized data set with mean and identity sample covariance. 

After standardization, one can now compare the standardized /-vectors to the standard 
Gaussian distribution, in terms of the Kolmogorov distance. The Kolmogorov distance be- 
tween two random variables X,Y is defined to be dx '■= sup tgK \ P[X < t] — P[Y < t]\. 
Observe that since the standardized data are conjectured to be approximately Gaussian, 
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any linear combination of their components should also be close to Gaussian. This moti- 
vates the following numerical test. 

Fix M, and let {di\ 1 ^L 1 C S p_1 be M independent samples of uniformly chosen points on 
the sphere in W . For each a iy compute the Kolmogorov distance 

<? K :=sup\F Nti (t)-$(t)\ 
teR 

where F N ^(t) is the empirical cumulative distribution function of the data projected onto 
the direction af, that is, 

and $(£) is the cumulative distribution function of the standard normal distribution. Finally 
define a measure D K of the distance from the data to Gaussian by 

D K := sup <f K . 

l<i<M 

The values presented in the tables in Section 1 are obtained using M = 100,000. 
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